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1.  Introduction 


The  feasibility  of  estimating  temperature  profiles  (synthetic  temperature 
profiles)  from  Geodetic  Earth  Orbiting  Satellite  (GEOSAT)  altimeter-derived 
sea-surface  heights  in  the  Gulf  Stream  region  has  been  explored  by  Carnes 
et  al  (1990).  The  scheme  was  based  on  a  statistical  relationship  between 
sea-surface  heights  (dynamic  height  at  the  surface  relative  to  1000  dbar)  and 
subsurface  temperature  profiles  derived  by  deWitt  (1987).  By  analysis  of  the 
U.S.  Navy’s  Master  Oceanographic  Observation  Data  Set  (MOODS)  hydrocast 
data  in  the  Gulf  Stream  and  Kuroshio  regions,  deWitt  found  that  the  first  two 
empirical  orthogonal  functions  (EOFs)  of  the  temperature  profiles  represented 
more  than  95  percent  of  the  overall  temperature  variance.  Furthermore,  he 
found  that  there  is  a  tight  relationship  between  relative  dynamic  height  and 
amplitude  of  the  first  two  EOF  modes.  This  relationship  was  used  by  Cames 
et  al  to  generate  the  synthetic  temperature  profiles  from  GEOSAT  data.  The 
synthetic  temperature  profiles  compared  well  with  the  expendable 
bathythermograph  (XBT)  measurements. 

In  this  study  the  temperature  and  salinity  profiles  in  the  Gulf  of  Mexico 
were  collected  and  analyzed.  Then  the  deWitt  scheme  was  employed.  The 
feasibility  of  using  sea-surface  heights  to  estimate  temperature  profiles  (synthetic 
profiles)  in  the  Gulf  of  Mexico  was  investigated. 

2 .  Data 

The  data  used  in  this  study  were  obtained  from  1)  a  quality-controlled 
version  of  the  National  Oceanographic  Data  Center  (NODC)  Oceanographic 
Station  Data  (SD)  file  (Levitus,  1982);  2)  the  U.S.  Navy’s  MOODS;  and 
3)  conductivity-temperature-depth  (CTD)  profiles  from  NODC  which  are  not  in 
the  SD  file  or  MOODS.  The  NODC  SD  file  contained  values  of  temperature, 
salinity,  and  dissolved  oxygen  at  NODC  standard  depths.  Temperature  and 
salinity  data  were  used  as  in  the  analysis. 

The  MOODS  contains  mostly  XBT  profiles  which  cannot  be  used  due 
to  1)  no  salinity  measurement  and  2)  no  data  deeper  than  700  meters.  Some 
CTD  profiles  were  available  in  the  MOODS  and  were  used  in  the  analysis. 
CTD  data  from  MOODS  were  quality-controlled  and  re-sampled  at 
NODC  standard  depths. 
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Some  CTD  profiles  which  were  not  in  the  SD  file  and  MOODS  were  also 
included  in  the  analysis.  Those  CTD  data  are  readily  quality-controlled. 
No  further  quality  control  was  done  except  a  re-sample  to  the 
NODC  standard  depths. 

A  subset  of  data  was  selected  for  the  region  of  the  Gulf  of  Mexico.  The 
hydrocast  profiles  that  were  less  than  1000  meter  depth  and/or  have  less  than 
10  valid  temperature  and  salinity  values  were  excluded.  A  total  of 
2388  hydrocast  profiles  were  available  for  the  analysis.  Figure  1  shows  the 
region  of  study  and  the  spatial  distribution.  The  data  were  divided  into  monthly 
subsets.  The  number  of  hydrocasts  used  at  each  standard  depth  is  listed  by  month 
in  Table  1.  The  monthly  mean  with  root-mean-square  (RMS)  deviation  are 
shown  in  Figure  2  for  temperature  and  in  Figure  3  for  salinity,  respectively. 

For  the  following  analysis,  the  temperature  and  salinity  at  19  depths  (0,  10, 
20,  30,  50,  75,  100,  125,  150,  200,  250,  300,  400,  500,  600,  700,  800,  900,  and 
1000  meters)  were  actually  used. 

3.  Empirical  Orthogonal  Decomposition 

From  each  monthly  temperature  data  set,  the  covariances,  Ry  where  i  and  j 
denoted  the  depths,  were  computed, 

Ry«£  I  (Tilt  -  Ti)  (Tjk  -  Tj)  i  =  1.  N;  j  =  1,  N, 

K  k=1  (1) 

where  N  is  the  number  of  levels,  K  is  the  number  of  hydrocasts,  and  T  is  the 
mean  temperature  profile.  The  covariance  matrix,  {Ry } ,  is  real  and  symmetric. 
Thus,  real  eigenvectors,  <t>n,  and  positive  eignvalues,  Xn>  can  be  found,  such  that 

N 

RijOni  —  A.nOnj;  n  =  1,  2,  •  •  -,  N, 

i  =  1  (4) 

where  the  number  of  modes  equals  the  number  of  depths.  The  resulting 
eigenvectors  are  orthogonal  and  are  referred  to  as  EOFs.  The  EOFs  of  the 
first  three  modes  are  shown  in  Figure  4. 
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Figure  1.  Region  of  Study  and  Spatial  Distribution  of  Hydrocast  Profiles 
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Table  1.  Number  of  Hydrographic  Observation 
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Figure  4a.  The  Empirical  Orthogonal  Function  (EOF)  of  the  First  Mode.  The  EOFs  represent 

decoupled  vertical  structure  of  temperature  variations. 
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Mgure  4b.  The  Empirical  Orthogonal  Function  (EOF 
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figure  4c.  The  Empirical  Orthogonal  Function  (EOF)  of  the  Third  Mode. 


The  contribution  of  each  mode  as  the  proportion  of  overall  variance  is 
determined  by  the  eigenvalues.  Table  2  lists  the  accumulate  contribution  in 
percentage.  The  first  two  modes  account  for  more  than  95  percent  of  the  overall 
variance  in  most  of  the  12  months. 

A  detailed  analysis  of  the  fit  by  reconstructing  the  temperature  profile 
using  various  number  of  modes,  such  that 

M 

Tj  =  Tj+  X  AnOnj;  n  =  1,  2,  •  • -M  and  M  <  N, 

n=l  (3) 


is  presented  in  Table  3.  The  misfit  represented  by  the  root-mean-square 
errors  (RMSE)  is  the  limit  of  accuracy  which  can  be  achieved  by  the 
EOFs  accounted  for. 


4.  Statistical  Relationship  between  Dynamic  Height  and  Modal 
Amplitudes 


The  relative  dynamic  heights  (0:1000  decibar)  were  computed  according  to 
the  hydrostatic  equation. 


(4) 


where  p  is  pressure,  p  is  density,  and  g  is  acceleration  due  to  gravity.  The 
amplitudes,  An,  of  the  first  three  EOF  modes  were  also  computed  from  cast  data. 


An  =  I  (Ti  -  Ti)  <t>ni; 


i  =  1 


for  n  =  1,  2,  3. 


(5) 


A  statistical  relationship  between  dynamic  height  at  surface,  T),  and  modal 
amplitudes  was  constructed  by  fitting  a  third-order  polynomial  regression 
equation, 

An  =  bon  +  binTl  +  b2n02  +  b3n  (6) 
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Table  3.  Analysis  of  the  Fit  by  Reconstructing  the  Temperature  Profile  Using  Various  Number 
of  Modes.  The  misfit  represented  by  the  root-mean-square  errors  (RMSF)  is  the  limit  of 
accuracy  which  can  be  achieved  by  the  EOFs  accounted  (November  shown  only). 


to  the  data  using  least  squares.  The  estimated  regressions  for  the  first  two 
EOF  modes  were  shown  in  Figure  5  for  November.  The  goodness-of-fit 
represented  by  R-squared  (in  percent)  is  very  high  for  the  first  mode  in  every 
month  (listed  in  Table  4).  For  the  second  mode,  the  goodness-of-fit  was  much 
lower  and  varied  each  month. 

5.  Synthetic  Temperature  Profiles 

Based  on  the  derived  statistical  relationship  between  dynamic  height  at 
surface  and  modal  amplitudes,  the  synthetic  temperature  profiles  were  estimated 
using  simulated  sea-surface  heights  and  GEOSAT-measured  sea-surface  heights. 
The  simulated  sea-surface  heights  were  computed  from  hydrocast  data  as  the 
relative  dynamic  height  (0:1000  decibar).  Results  of  error  analysis  of  the 
synthetic  temperature  profiles  from  simulated  sea-surface  heights  are  listed  in 
Table  5.  One  section  of  synthetic  temperature  profiles  along  the  GEOSAT  track 
(Figure  6)  is  demonstrated  in  Figure  7. 

6.  Summary 

The  statistical  relationship  between  modal  amplitudes  and  surface  dynamic 
height  in  the  Gulf  of  Mexico  has  been  investigated.  Monthly  temperature  data 
were  used  to  determine  the  EOFs  and  the  modal  amplitudes.  Temperature  and 
salinity  data  were  applied  to  computed  sea-surface  heights  (relative  dynamic 
heights).  Synthetic  temperature  profiles  were  generated  using  simulated 
sea-surface  heights  and  the  related  errors  were  estimated.  Synthetic  temperature 
profiles  generated  from  GEOSAT-measured  sea-surface  heights  were  also 
demonstrated. 
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Figure  5.  Statistical  Relationship  Hetwcen  Modal  Amplitudes  of  KOFs  and  Dynamic  Height. 
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Table  4.  The  Regression  Coefficicncies  R-Squarcd  (in  percent)  for  the  First  Two  Modes. 
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Table  5.  The  Error  Analysis  of  the  Synthetic  Temperature  Profiles  From  Simulated  Sea-Surface 

Heights  (November  shown  only). 
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Figure  6.  GEOSAT  Exact  Repeat  Mission  (ERM)  tracks  over  the  Gulf  of  Mexico. 
Along  track  sea-surface  height  variations  (solid  line)  were 
used  to  generate  synthetic  temperature  profiles. 
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K  Synthetic  Temperature  Profiles  on  a  Cross  Section  Along 
GEOS  AT  ERM  Track  (Shown  by  solid  line,  Figure  6). 
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